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ABSTRACT 



Ar experimental and theoretical analysis of a rotating system with 
flexibly supported, viscously damped bearings is described. The theo- 
retical development is based on the differential equations of motion of 
the mass of the rotor and the two shaft bearings. Tests were conducted 
and compared with calculated response of the system. Results indicate 
fair agreement The discrepancies are attributed to limited accuracy 
of test measurements and/or erratic operation of the viscous damping 
dashpo ts o 

This work was conducted at the U« S Naval Postgraduate School, 
Monterey, California . 
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1 INTRODUCTION 



The problem of vibration in rotating systems has commanded the 
attention of theoreticians and practical investigators for many years. 
With the advent of high speed marine turbines and turbo-jet aircraft 
engines renewed interest has been shown in this field. 

It is well known that at certain speeds of rotation small imbalances 
tend to set up vibrations of. large amplitudes in a rotating system, 

Prohl, who has made several contributions in this area, ^ 1 , 2 | , wrote, 

The existing methods for determining critical speeds are 
subject to the following limitations. On one hand the methods 
that are general, i e. , that permit the calculations of higher 
speeds as well as the fundamental, involve calculations so 
complicated as to be impractical for any but the simplest of 
rotors. On the other hand, the methods for which the calculations 
are compare t ively simple, such as the more familiar methods of 
Rayleigh and Stodola, lack generality in that critical speeds 
other than the fundamental cannot be definitely determined. 

This criticism may have been well taken at the time. The large 
number of prime variables in any but the simplest system tend to confine 
most investigations to a very narrow scope; for example, M 0n Critical 
Speeds of a Shaft Supported by a Ball Bearing 11 , by T. Yamamoto, which 
appeared in the Journal of Applied Mechanics of June, 1959. Doctor 
Yamamoto investigated two kinds of critical speeds induced by a slight 
difference in the diameter of the balls in their races. 

A growing appreciation of modern computer capabilities has led to 
the development of increasingly general treatments such as, M Natural 
Frequencies of Nonuniform Beams on Multiple Elastic Supports 11 , by 
R A DiTaranto of. the Defense Electronics Products Division of Radio 
Corporation of America, | 3j. An especially general and pertinent 
method was described recently by Messrs „ Koenig, Guenther and Love joy 
of the Engineering Analysis Section, Ail is -Chalmers Manufacturing 
Company, \ 4j A concluding, third part of this method, as yet unpub- 
lished, is to include experimental results and correlation with the 
analytical method described 
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The study presented here investigates the transverse motion of 
a rotor on a flexible shaft- The motion of the two flexibly supported 
end bearings, which are subject to damping, is also studied. 

The study includes both a practical investigation conducted with 
a machine designed specifically for this purpose and a theoretical 
analysis developed along parallel lines. 

The theoretical analysis begins with the basic equations of motion. 
This approach was considered preferable to a mechanical application of 
one of the existing schemes of analysis. 

The analytical model was constructed with the assumption that the 
rotor was the only rotating member. The shaft was considered to be 
massless and torsionally rigid The speed of rotation was taken to 
be constant. Damping was assumed to be a linear function of the 
velocity 

This work was conducted at the U. S. Naval Postgraduate School, 
Monterey, California during the 1961-1962 academic year. 

The author wishes to express his appreciation to Professors 
T. H. Gawain and H. Vavra and to the supervisory and laboratory 
personnel of the Department of Aeronautics of the U. S. Naval Post- 
graduate School. 
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2 TABLE OF SYMBOLS 



A general constant 

Fortran array 

a rotor position measured from left-hand support, in. 

C (1, 2, etc,) Fortran physical constant 

(A, E) Fortran damping constant 

c symmetrical damping constant, lbs /(in. /sec.) 

D displacement, in 

mean spring diameter, in. 

d wire diameter, in. 

F force, lbs, 

G shear modulus, lbs./sq, in, 

2 

g gravitational constant, 386 in.. /sec. 

4 

1 diametral moment of inertia, in 

i imaginary, (~1) 2 

K,k spring constant, lbs./in, 

X shaft length, in, 

2 

m mass, lbs. - sec, /in, 

N number of working coils of spring 

first critical speed, RPM 

P load, lbs . 

Q general constant 

R general constant 

t time, sec, 

W weight, lbs* 

w rim weight lbs 

x displacement of rotor c.g, from polar axis, in. 

x displacement, in. 

x velocity, in. /sec 

, . . , 2 

x acceleration, m /sec 

Z general constant 
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Table of Symbols 



6 


complex root 


e 


rotor eccentricity, in. 


\ 


damping factor, sec. ^ 




mass ratio 


*T 


“l + P, 3 


5 


*1 ^3 


tt 


3,1416 


0 


phase angle, degrees 


(JU 


angular frequency, radians/ 


Subscripts 




1 


left hand bearing complex 


2 ,R 


ro tor 


3 


right-hand bearing complex 


d 


damping 


I 


inertia 


T 


total 


Abbreviations : 


c.g 


center of gravity 


cps 


cycles per second 


in. 


inches 


lbs . 


pounds 


m , u o 


mass units 


RPM 


revolutions per minute 


sec . 


seconds 


sq . in . 


square inches 
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3. THE PHYSICAL SYSTEM 



A schematic drawing of the physical system is shown in Fig. 1. 
Detailed drawings from which the critical speed test rig was constructed 
are included in Appendix A. 

3 1 Critical Speed R i g Design Concept 

Bearings and restraints of high-speed rotating systems should be 
designed with great care. Bearings are normally confined as rigidly 
as possible in order to reduce the degrees of freedom of the system and 
to eliminate critical speeds below the operating speed range. However, 
shaft and bearing flexibility are inevitable due to the elastic nature 
of materials, oil film "s ti f fnesses 1 ' , necessary clearances and machining 
tolerances. Similarly, imbalances due to imperfections and lack of 
homogeneity of shaft and rotor materials are always present to some 
degree. Even after exhaustive balancing procedures have been carried out 
some imbalances remain and contribute harmonic excitations in rotating 
sys terns . 

The term impedance is sometimes used in describing the character- 
istics of spring and damping device combinations which act at the supports 
of a system. As noted by Guenther and Lovejoy, [4], support mechanical 
impedance can be, and is often, expressed as a function of not only the 
static properties of a system but of the operating speed as well. 

In the case where support is provided by a simple spring, support 
stiffness is merely the familiar spring constant (force per unit 
deflection). The ideal mechanical spring acts as a conservative member 
in the system; energy stored in compression or tension is returned to 
the system when the spring returns to its unstrained state. Hysteresis 
losses present in actual springs are generally small and can be disre- 
gardedo 

For the simple ? damped, spring -mass system as discussed in most 
basic vibration texts, e.g,, [ 5], a dashpot is shown as a device which 
develops forces opposing motion, the force so developed being propor- 
tional to the velocity. A dashpot acts as a non-conservative member 
dissipating energy from the system. 
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Both the spring and the dashpot influence the critical speed(s) . 
The dashpot governs the amplitudes associated with critical speed 
vibra tions 



The critical speed test rig was designed to provide a vehicle for 
practical investigation of the effects of support flexibility and 
viscous damping on the critical speeds and mass vibration amplitudes of 
a rotating system. 

The salient points of the design with respect to its particular 
service are: (a) compressed air drive as opposed to mechanical drive 
which would constrain the shaft during vibration, and, (b) ability to 
control dashpot stiffness by means of a suitable valve. 

Choice in arrangement of the three support frames, the several 
rotor and spring sizes available, range in dashpot stiffness and choice 
of rotor position on the shaft afford a broad selection of configurations 
possible of investigation with this test rig. 

3. 2 Test Configuration and Major Components 

3*2.1 Drive Disk and the Bearing and Linkage Complex 
Drawings of these components are shown in Appendix A. Pertinent 
properties are listed in Table I, Physical Properties of Components. 
Figure 2 shows the drive assembly and Fig. 3, the bearing and linkage 
complex* 



3.2.2 Springs 

Figure 3 shows springs mounted in support of one of the shaft 
bearings. Details of the various springs are shown in Appendix A and 
listed in Table 1. 



Spring constants weie calculated by formula, i 6 
P = load, lbs. 

F = deflection caused by given load, in. 

G = shear modulus, lbs^/sq. in. 
d = wire diameter, in, 

D = mean diameter of spring, in. 

N = number of effective working coils 



P 

F 



G d 



8 N D 



3 



lbs /in. 
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PHYSICAL PROPERTIES OF COMPONENTS 
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Cusped, aluminum drive disk 
Compressed air nozzle 
Control v-°lve 




Bearing ring 
Linkage to drshpot 



Light Spring 



P = 1 'i , 4 x 1 0 & x 2 0 7 x 10 

F 8x12x1.4 



17.5 lbs. /in. 



Intermediate Spring 

P = 11.4 x IQ 6 x 9.82 x 10 

F 8 x 8 x 1 . 63 



107 lbs. /in. 



Spring constants were also determined experimentally. The spring 
and its bolder were set vertically over a slot in a fixed surface. A 
small wire attached to the top holder was led along the axis of the 
spring and through the slot in the table. The unloaded length of the 
spring was measured with a vernier caliper. Weights of known sizes 
were attached to the free end of the wire. The length of the spring 
was measured for each increment of step-wise loading, and similarly 
during unloading Results are listed in Table II, Calibration of Light 
and Intermediate Springs. Results are shown graphically in Figs. 4 and 5 
for light and intermediate spring respectively. Calculated and measured 
values of spring constants are listed in Table III, Measured and Calcu- 
lated Spring Constants, for purposes of comparison. 

The lack of agreement in the calculated and test values for the 
light spring indicates that friction in the holders may have had a strong 
influence on the experimental results. Further testing of the spring is 
advisable Both Fig. 4 and 5 show pronounced hysteresis losses are 
present during compression and relaxation of the springs. The slopes of 
the figures, from which the response was calculated, were not strongly 
effected by hysteresis. The experimental value of the intermediate spring 
constant, in good agreement with the calculated value, was used in all 
subsequent calculations Only the intermediate spring was used in tests 
for critical speeds • 
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Table II 



w> 



CALIBRATION OF LIGHT AND INTERMEDIATE SPRINGS 





LIGHT SPRING 


INTERMEDIATE SPRINGS 


Wt. (lbs. 


) Def lection(in. ) 


Wt. (lbs. ) 


Def lection (In . ) 


0.0 


0.0 


0.0 


0.0 


1.0 


0.053 


1.0 


0.061 


2.0 


0.062 


2.0 


0.109 


3.0 


0.066 


3.0 


0.152 


4.0 


0.078 


4.0 


0.238 


5.0 


0.083 


5.0 


0.290 


6.0 


0.091 


6.0 


0.382 


7.0 


0. 121 


7.0 


0.513 


8.0 


0.135 


8.0 


0.597 


9.0 


0,146 


7.0 


0.550 


10.0 


0.168 


6.0 


0.506 


9.0 


0.149 


5.0 


0.447 


SO 


0-143 


4.0 


0.365 


7.0 


0.133 


3.0 


0.281 


6.0 


0. 122 


2.0 


0.214 


5.0 


0.110 


1.0 


0.142 


4.0 


0. 101 


0.0 


0.101 


3.0 


0.087 






2.0 


0.082 






1.0 


0.070 






0.0 


0.062 
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Table III 

MEASURED AND CALCULATED SPRING CONSTANTS 



Spring 


Measured 

Constant 


Calculated 

Constant 


7. Variation 


Light 


14.5 


17.5 


20. 5 


Intermediate 


105.5 


107.0 


1.4 






* 

based on 


measured value 



3.2,3 Shaft 

Dimensions and physical properties are shown in Appendix A and 
Table I. 

The shaft was checked for set or eccentricity by rolling on a 
precision flat surface. Variations noted were so slight that they did 
not warrant correction. 

Static shaft spring properties were calculated for midspan loading 
with pin-ended supports. The length of the shaft between supports was 

3 

taken as 40*5 inches. From elementary beam theory, K = 48 El/ (length) 
lbs. /in. The calculated value was 65.33 lbs. /in. 

The spring properties of the shaft were also evaluated experimentally. 
Tests were made with the shaft on knife-edged end supports and with the 
shaft supported in the bearings. Any change in response would indicate 
a resisting moment could be exerted despite the self-aligning feature of 
the bearings. Tests were made for the unloaded shaft and with the rotor 
in position at midspan. With the shaft on supports placed 40.5 inches 
apart loading and unloading was conducted in the same manner described 
for the springs. Results are listed in Tables IV and V for the knife- 
edge and bearing support end conditions respectively. Figure 6 shows 
the experimental results graphically. No significant change resulted 
from the change in end support conditions. From the slope of the curve 
in Fig. 6 the shaft stiffness was calculated as 66.5 lbs. /in. Although 
the test with the bearing showed some hysteresis loss the effect was 
small, and was therefore disregarded. The experimental value, in close 
agreement with the theoretical value, was used in all subsequent cal- 
culations c 
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Table IV 



STATIC SHAFT SPRING CALIBRATION - KNIFE EDGE SUPPORTS 



Wt. (lbs.) 


Deflect ion(in. ) 


*W „ 
rotor 

Wt . ( lbs . ) 


Def lection(in. ) 


0.0 


0.0 


W +0.0 


0.0 


1.0 


0.014 


r +1.0 


0.016 


2.0 


0.029 


+2.0 


0.032 


3.0 


0.042 


+3.0 


0.046 


4.0 


0.057 


+4.0 


0.064 


5.0 


0.072 


+5.0 


0.080 


6.0 


0.088 


+6.0 


0.090 


7.0 


0.103 


+7.0 


0.107 


8.0 


0.116 


+6.0 


0.093 


9.0 


0.132 


+5.0 


0.079 


10.0 


0.146 


+4.0 


0.064 


11.0 


0.160 


+3.0 


0.048 


12.0 


0.178 


+2.0 


0.032 


11.0 


0. 164 


+1.0 


0.016 


10.0 


0.147 


W +0.0 


0.003 


9.0 


0. 132 


r 




8.0 


0.117 






7.0 


0.102 






6.0 


0.086 






5.0 


0.072 






4.0 


0.056 






3.0 


0.042 






2.0 


0.029 






1.0 


0.014 






0.0 


0.0 







*W = 14.5 lbs. 
r 
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Table V 



STATIC SHAFT SPRING CALIBRATION - BEARING SUPPORTS 



Wt (lbs.) 


Def lection(in. ) 


*W 

rotor 
Wt. (lbs. ) 


Deflection(in. ) 


0.0 


0.0 


W +0.0 


0.0 


1.0 


0.015 


r +1.0 


0.019 


2.0 


0.028 


+2.0 


0.041 


3.0 


0.041 


+3.0 


0.050 


4.0 


0.061 


+4.0 


0.070 


5.0 


0.076 


+5.0 


0.084 


6.0 


0.091 


+6.0 


0. 103 


7.0 


0.108 


+7.0 


0.115 


8.0 


0.125 


+6.0 


0.103 


9.0 


0.145 


+5.0 


0.090 


10.0 


0.157 


+4.0 


0.074 


11.0 


0.172 


+3.0 


0.060 


12.0 


0.190 


+2.0 


0.046 


11.0 


0.181 


+1.0 


0.030 


10.0 


0.164 


W +0.0 


0.014 


9.0 


0.149 


r 




8.0 


0. 135 






7.0 


0.121 






6.0 


0. 103 






5.0 


0.088 






4.0 


0.076 






3.0 


0.055 






2.0 


0.042 






1.0 


0.029 






0.0 


0.016 







*W =14.5 lbs. 
r 
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3.2.4 Rotor and Lip Weight 

Details of the rotor are shown in Appendix A and listed in Table 1. 
Figure 7 shows the relative size of these parts. 

The rotor was balanced statically. Dynamic balancing was considered, 
but it was felt that the disk, only 1.5 inches thick, would produce 
only small effects due to possible variations in properties between its 
parallel planes. 

The rim weight, which was introduced after static balancing, was 
employed for two purposes. Knowledge of its size permitted calculation 
of rotor eccentricity, defined as the shift in center of gravity of the 
rotor due to the. added weight. This value is essential to accurate 
calculation of exciting forces developed during rotation. Secondly, 
considering the aforementioned static balancing, the imbalance due to 
the lip weight would overshadow effects of any small variations not 
removed in the balancing. 

The eccentricity was calculated as follows: 

W = basic rotor weight, lbs. 

w = added rim weight, 0 04 lbs @ 2.75 in. from axis 
x = location of displaced c.g. from polar axis, in. 

Moment about x = W x ~ w(2.75 ~ x) = 0 



0.04 x 2.75 
14.44 + 0.04 



0.0076 in. 



e 



3*2.5 Control System 

Details of the components of the control system are shown in 
Appendix A The general arrangement is shown in Fig. 2. 

The item referred to as the control valve in Fig. 2 was not designed 
specifically for this use. It is identical to the needle valves used as 
metering devices for the working fluid in the dashpots. Large adjustments 
in supply of compressed air to the drive jet were made by means of a 
coarse valve in the line upstream of the control valve. 

32,6 Dashpots 

Detailed drawings of the dashpots are shown in Appendix A. 

Smith, ^7], has written, 

The assumption as to the viscous law of damping is not fulfilled 

in practice and analytical predictions as to the effect of damping 

are, to some degree, only qualitative. 
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Figure 7 Rotor ona Lip /eight 



The dashpots required extensive testing to determine their char- 
acteristics, It was desired to determine three main factors: 

(1) the degree to which the dashpots were linear, (2) the dependence 
of the magnitude of damping on the direct communication port openings 
(See Fig. 8), and (3) the range of damping "constants 11 possible 
within limits of needle valve settings. 

Any attempt to evaluate the dashpots by analysis of the internal 
flow would require numerous assumptions. An analytical description of 
the non-steady flows through the bypass tubing and direct communication 
ports would have resulted in a doubtful approximation at best. 

A method based on direct measurement of drive force, acceleration 
and velocity response was used. 

3.2. 6.1 Instrumentation and Test of Dashpots 

A dashpot was mounted on a fixed platform above a Ling Electronics 
Model 219 Shaker as shown in Fig. 9. A drive spindle of 6061-T6 aluminum, 
\ 8 ], was made to accept the male thread of the dashpot piston shaft. The 
spindle, except for the end portions, was square in cross section. The 
width of the sides of the cross section was machined to 0.250 + 0.001 
inches. Two SR-4, type A- 7 strain gages, (9], were mounted axially on 
opposite sides of the spindle. The gages were connected in series to 
increase sensitivity to strain, and at the midpoint of the spindle to 
avoid effects from the non-uniform ends. The other end of the spindle 
was threaded into an adaptor on the drive table of the shaker. 

Two accelerometers were mounted on the adaptor. Figures 10 and 11 
illustrate the details of the arrangement. A Statham Model 5A5-100-120 
accelerometer was used to measure accelerations. Its output, as well 
as that from the strain gages, was led to the Honeywell Model 130-2C 
Carrier Amplifier shown as item 7 in Fig, 12. Both amplified signals 
were recorded semi.” permanently on the output paper of the Heiland 
Model 906B Oscillograph Visicorder, item 6, Fig. 12. An electronic 
timer imposed time reference lines on the simultaneous trace of accel- 
eration and strain. The output of the second accelerometer, a Consol- 
idated Engineering Corporation Model 4-102-A Vibration Pickup, measured 
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Figure 8 Doshpot Piston 
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velocity and peak* to -peak amplitude of vibrations. Readings of average 
velocity and displacement were indicated on a Consolidated Electro- 
dynamics Corporation Type 1-117 Vibration meter. The velocity, amplified 
by the vibration meter, was led to a Hewlitt-Packard Model 120A 
Oscilloscope A typical trace of the variation of velocity with time 
is shown in Fig. 13. Operating frequency was read directly from the 
frequency control dial on the shaker control console. 

The strain trace on the Visicorder was calibrated in accordance 
with operating instructions for the system. The dashpot was charged 
with H-.515 OHA petroleum base hydraulic fluid, [lOj. Several combin- 
ations of port openings and bypass valve settings were observed. 

Operation with hydraulic fluid in the dashpots was conducted in some 
cases and air was used in other instances. 

It. was believed that even the minimum possible damping on the part 
of the dashpots might prove excessive in relation to other characteristics 
of the system, It was decided to attempt to approximate the case of 
zero damping. 

The dashpot configuration which would produce minimum damping 
obviously required the bypass line and communication ports to be unre- 
stricted. The only question lay in the choice between air and hydraulic 
oil as the working fluid. Although the choice of air would seem the 
more logical, there was a higher likelihood of large frictional forces 
associated with its use. It is not known whether an air loaded dashpot 
would, in fact, produce less damping than that which was measured, but 
the hydraulic fluid was selected for use in the test based on its lubri- 
cating properties . 

The plugs were removed from the communication ports and the bypass 
valve was set fully open. Average velocity for the test was 3.0 in. /sec. 
Peak- to -peak displacement was 0c28 in. at the test frequency of 5 cps. 

Traces of acceleration and strain (force) were scaled up four-to-one 
for readability. Figure 14 shows only the basic variations of force 
and accelerations as recorded on the trace. Superimposed on the actual 
trace was a high frequency M hash n . This was present during the warm up 
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Figure 13 Piston Velocity Trace 
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of the shaker prior to actual tests. Since these signals were obviously 
not associated with response of the dashpot they were disregarded in the 
interpretation of the test traces. By successive integrations of the 
acceleration trace, velocity and displacement traces were produced. 

These waveforms are shown in Fig. 14 a-c. The associated force trace 
is shown in Fig^ 14 d. Integration of the acceleration wave produced a 
parabolic velocity trace. The velocity trace was scaled from the known 
value of the average velocity, two- thirds of the maximum velocity repre- 
senting the average height of the parabolic trace. The scale of the 
acceleration trace was found from the fact that the area under the accel- 
eration curve was equal to the known maximum velocity. 

According to design specifications the shaker should have produced 
a sinusoidal output with maximum of 3 percent distortion at 5 cps . The 
derived displacement- time relationship was not a sine wave. Whether this 
was due to instrumentation error or to shaker operation with distortion 
greater than that cited in the specifications is not known. 

With the assumption of sinusoidal motion for the system depicted in Fig. 15 
the following relationships apply: 

Fj = inertia force = m x ; F<^ * damping force = c x 

F T = total force = F^. + F^ 

2 

x = A sin UDt , x - cx)A cos out , x = - (i) A sin (Dt 






F^ * c x , a straight line whose slope is c. 

2 2 
Fj = - m U) x , a straight line whose slope is - m to . 
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Ideal Graphic Relationships 



ctual 



Actual 




Figure 15 shows these relationships graphically. The ellipse in the 
velocity plot has a minor semi-axis proportional to the inertia or 
mass of the system The minor semi-axis of the displacement plot is 
proportional to the force of viscous damping in the system. Repre- 
sentative plots of experimental data are included for purposes of 
comparison 

Consistent with the nature of the forces involved the inertia force 
contributes nothing to the area of the displacement plot. 

Test data plotted in Fig. 16 indicate the variation of total force 
with velocity. The shape of the enclosed area differs markedly from the 
ideal form of an ellipse. The distortion of the figure emphasizes the 
remarks cited at the beginning of this section. The non-linearity is 
strongest in the region of negative velocity. Deviation from the ideal 
case is also apparent in the recorded force trace. Fig. 14 d. This may 
have been caused by coverage of the bypass ports at the extreme point of 
the piston stroke, by large friction forces present due to momentary 
loss of lubrication between the dashpot cylinder and the piston surfaces, 
or by lag in the mechanical linkage between the drive table and the 
piston. The vertical offset of the center of the closed figure from the 
origin indicates that the force is dependent to some degree on the 
absolute direction of motion. No logical reason for such a dependence 
is apparent. 

For a linear system the slope of the mean line of the F^ - velocity 
plot is equal to the damping constant. Since no dependable mean line 
is apparent, the slope of a line drawn through the endpoints of the figure 
is considered only as a first estimate. The slope so defined was cal- 
culated as 1.88 lbs . /(in. /sec. ) 

In an attempt to analyze the F^ - velocity trace further, the product 
of mass and acceleration was subtracted point-by-point from the total 
force ordinate. Figure 17 shows the result as a "Reduced Force"-velocity 
plotc A further refinement of this figure could be made by an allowance 
for the mass of fluid being accelerated during the operation of the 
dashpot, Regardless of the distortions in Figs. 16 and 17 it can be said 
that, essentially, the dashpot response has been established. 
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The force-displacement plot from test data is shown in Fig. 18. 

The closed curve has some resemblance to an ellipse. Assuming that 
an equivalent damping system will absorb or dissipate equal energy per 
cycle as does the physical system, an equal area ellipse was constructed. 

One semi-axis was defined by the known displacement. From the formula 
for the area of an ellipse, Area = tt x a x b, the regular ellipse, shown 
dashed in the figure, was drawn and superimposed. From the scale of the 
test figure the semi-axis proportional to c(DA was measured. Dividing by 
the known value of (DA the value of the damping constant was found. 

c ’ —vHVx'o.'w ' 1 ' 77 H>»-/(in./s.c> 

This value, in fair agreement with that calculated from Fig. 16, is 
probably the better of the two. 

3. 3 Critical Speed Test Program 

Tests of the rotating system were conducted with the bearings locked 
and with the bearings under the restraint of the dashpot-spring combination. 
The dashpots for the latter test were set in the same configuration used 
in their calibration, i.e., direct communication ports and bypass control 
valve fully open. The intermediate spring was employed with the dashpot. 

3.3.1 Instrumentation 

The electrical probe displacement sensor shown in Fig. 19 was used 
to measure deflections. A conducting probe in the form of a small 
spring was fixed to the end of the vernier shank. The vernier micrometer 
was mounted rigidly in a support. A 22.5 volt dry cell battery provided 
potential to the probe through the conductor shown. This voltage was 
impressed in parallel across the vertical deflection circuit of the 120A 
Oscilloscope. At each instant of contact between the probe and the 
rotor the scope circuit was grounded and a vertical spike appeared on the 
scope. Zero readings of the vernier were made prior to each test. The 
vernier, readable to 0.001 inches, was then backed off. For each speed 
at which measurement of vibration amplitude was desired the vernier was 
advanced until first contact was indicated on the scope. Some arcing 
at the tip of the probe indicated small errors, estimated at + 0.005 in., 
were inherent to the method. 



- 35 - 






3± 


+j 


fj 




if 




-rf 


H 


-fr 


3 


5 


p 


Tf 

TT 


-rr" 


-LJ. 


4-f 


±t 


# 










Figure 19 Electric Probe Displacement Sensor 




An RCA Model 631B Strobotac was used for speed indications. The 
low scale on this instrument, usable for speeds from 60 to 3700 RPM, was 
readable to + 2 RPM. 

3.3.2 Results 

At no time during tests with the bearings supported by the inter- 
mediate spring was there any perceptible or measurable motion of the 
bea rings • 

Results are listed in Table VI, Experimental Results. Results are 
shown graphically in Figs. 31 and 39 in conjunction with the comparison 
of theoretical and experimental results. 

3.4 Comments on the Physical System 

The system could not be tested without the dashpot linkage, and 
consequently, the d&shpots engaged. This resulted from the inability of 
the spring-spring holder combination to resist loads normal to its 
axis. The stabilization provided by the dashpot linkages was desirable 
in one sense and indicative of a possible source of trouble in another. 
Forces along the rotor shaft axis were carried through the dashpot 
linkages to the dashpot pistons. An aggravated case of this could have 
resulted in the pistons becoming cocked in the dashpot cylinders. Such 
cocking could have caused deviations from linearity even greater than 
those observed in the controlled testing of the dashpots. 

The light spring was not capable of supporting the combined weight 
of its associated bearing and the rotor mounted on the shaft. These 
springs "bottomed" in their holders when the fourteen pound rotor was 
mounted. A lighter rotor, capable of being supported by the light 
springs, was mounted. Spring action was linear with the rotor at rest 
but, movement of the bearings in vibration closed the springs coil-to- 
co il . 

The drive disk whipped noticeably during strong vibration of the 
system. Air was spilled from the cusps in the disk and a non-steady 
drive resulted The disk could be reworked to form channels as in a 
Pelton type turbine; this would reduce spilling. This same modification 
would eliminate the small axial component associated with the present 
drive system. 
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260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

372 

375 

380 

385 

390 

395 

400 

450 

500 



Table VI 



EXPERIMENTAL RESULTS 



fixed Hydraulic fluid in dashpots, 

bypass and ports open 



Deflection (in.) 


RPM 


Deflection (in.) 


0.044 


300 


0.001 


0.045 


330 


0.000 


0.046 


340 


0.005 


0.047 


355 


0.029 


0.048 


370 


0.040 


0.050 


375 


0.080 


0. 062 


380 


0. 104 


0.063 


383 


0. 116 


0.066 


388 


0. 109 


0.067 


390 


0.113 


0.090 


393 


0.090 


0.099 


400 


0.075 


0.2 38 


450 


0.058 



0.670 
0.969 
1 . 245 
0. 982 
0.072 
0.051 
0 017 
0 004 
0 003 
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Strong transient vibrations were noted during accelerations of the 
rotor. The speed control system is not conducive to setting of precise 
speeds. Low speeds require small inputs from the jet. Acceleration forces 
must be kept very close to those necessary for steady drive at the speed 
desired or the speed will build rapidly and overshootc For low speeds 
this means allowance for long periods of speed buildup Fifteen minutes 
was not excessive for speed stabilization during tests. 

Fluctuations in the supply air pressure caused speed to drift even 
after long periods allowed for stabilization. In consideration of these 
factors measured speeds can be interpreted as accurate only to within 
+ 5 RPM. 

Error in the calculation of the dashpot damping constant may have 
stemmed from either or both of the following: (1) Low amplification of 

the spindle strain trace such that its maximum value reached only twenty 
percent of the lowest mark on the calibration scale of the recorder. 

(2) Reading and scaling of recorded traces. These sources of possible 
error could be eliminated or minimized by, (1) use of more sensitive 
strain gages with larger output per unit strain, and, (2) recording 
velocity waveform in addition to simultaneous acceleration and drive 
spindle strain. This would require the use of a strain gage type 
vibration pickup^ the Consolidated instrument used in the tests described 
was a self -exciting , piezo-electric type whose output could not be 
matched to the recorder galvanometer. 

Other configurations of bearing support, in particular that with the 
light spring mounted in the horizontal plane, were observed- Bearing 
motion was pronounced and the spring appeared to be operating without 
bottoming. This constituted a system with non- symmetrical bearing 
stiffness. Although not within the scope of this study, this configuration 
is of considerable interest any may form a logical extension to this work. 
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A . THEORETICAL ANALYSIS 



A model representative of the rotating system is shown in Fig. 20 
Component names and general nomenclature associated with the physical 
system were retained wherever possible. 

Values of the physical properties of the system used in this 
analysis differ somewhat from those listed in Table I. Only three 
masses are defined in the theoretical model; the mass of the rotor, 
m 2 , and the masses of the left and right-hand bearing complexes, m^ 
and m^, respectively 

A . 1 Equivalent Mass Lumping 

Any physical system can be said to have infinite degrees of freedom 
in that it is composed of an infinite number of point masses. In some 
cases less significant masses are ignored completely. An example is the 
elementary spring-mass system. The center of gravity of the spring is 
obviously displaced during oscillations and it therefore constitutes a 
mass in motion. Since the principal mass of the system is very often 
several orders of magnitude greater than that of the spring and the 
displacement of the spring mass is so small its contribution to the 
system is insignificant For the system under consideration here, an 
attempt is made to formulate a logical expression to reflect the con- 
tribution of the motion of the shaft to the total energy of the system. 
The kinetic energy of the actual shaft in vibration will be referred to 
as the shaft energy. The shaft energy will vary for each mode of 
vibration. Figure 21 shows, in a representative manner, the three modes 
of vibration for the three mass system. It is easily seen that no 
single factor will suffice in apportioning the mass of the shaft among 
the masses in the model. The first mode shown requires that the con- 
tribution be made solely to the rotor. The other modes indicate that 
the bearing masses should be increased by some factor but that the rotor 
mass need not be modified. If the elastic line of the shaft is assumed 
to be parabolic, eight-fifteenths of the shaft mass should be added to 
that of the rotor for similarity If the other modes were calculated 
on this basis they would be in error. Ideally an assumption should be 
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made as to the apportionment for each mode, the mode shapes calculated 
and the respective factors revised. Due to the relatively small mass 
associated with the actual shaft an iteration of this kind was not felt 
to be necessary. A factor of one-third was taken; the mass of each 
bearing and that of the rotor was increased by one- third the mass of the 
adjacent shaft section (s). The mass of the overhung drive disk, being so 
physically close to the left-hand bearing, was lumped with that bearing 
complex 

Table VII * Equivalent Lumped Masses, lists the calculated values 
of the three masses in the theoretical model. 

4 . 2 Mass Dynamics 

Free body diagrams of the three masses with forces acting during 
rotation are shown in Fig, 22. The differential equations of motion 
based on the free body diagrams reduce to a system of six simultaneous 
algebraic equations. The development is shown on the pages following 
Fig, 22. 

Had a solution in the form = ILe^ been assumed, the system 
could have been described in three equations with complex coefficients. 
Computer programming of complex expressions requires a somewhat more 
demanding technique, and, for this reason, the expressions with real 
quantities throughout were selected. 

The six equations describe the theoretical model of the rotating 
system in damped, forced vibration. As shown in the last page of the 
development, the system of general equations reduce to a specific case 
for the rotor at the middle of the shaft span. 

4.2.1 Program CRTSPD 

The six simultaneous algebraic equations are the basis of the 
computer program CRTSPD (Critical Speed). For tran 1 1 ! , coding was 
used. A program flow chart indicates the logic and progress of the 
computations in the computer. A sample program is included in 
Appendix B 
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Rotor mass: 



m R = + -y (.0544/386) (20 25) = 0.394 

Left-hand bearing complex mass: 

mi = ^ff + l ( 0544/386) (20.25) + 2 035 = 0 0270 

1 JC O J 

Right-hand bearing complex mass: 

m 3 = 7.88 + 1 (0544/386) (20.25) = 0.0220 



Table VII 

EQUIVALENT LUMPED MASSES 

m l 
m 2 

m 3 

Mass Ratios: 

m 1 /m 2 = 0.0270/0 394 = 0.6858 = (J. 
m 3 /m 2 = 0.0220/0.394 = 0.5588 = 



Significant Mass 
L. H. Bearing complex 

Ro tor 

R< H. Bearing complex 



W(lbs . ) 
10.434 

15.213 

8.501 



M(m . u . ) 
0.0270 

0.0394 

0.0220 
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Figure 23 



Flow Chart PROGR/i: CRTSPD 




The dimension statement apportions storage space for the arrays 

representing the augmented matrix of coefficients of the algebraic 

equations of motion and the values of their solution, X(6). Column 

headings for the printout are specified. The constants of the physical 

2 2 2 

system, C1-C6, represent oo , (JO , to , |JL , and e, respectively. 

R 1 J 1 j 

CA is the coded form of the damping constant at the left-hand bearing. 

CE is the corresponding damping constant for the right-hand bearing. 

These constants were set equal to one another for all calculations. The 
value so assigned is hereafter referred to as the symmetrical damping 
constant. The DO 10 statement defines the range of RPM. Omega is then 
defined in terms of RPM. Subroutine CALCA, to which all factors necessary 
for the calculation of the elements of the augmented matrix have been 
made available through use of the COMMON statement, then calculates each 
element according to the specified value of RPM. The resulting numerical 
array, A, common to Subroutine J0RDAN2 , is solved as a set of simultaneous 
equations. The values of the solution are returned to the main program 
and amplitudes of mass motion and associated phase angles are calculated. 

Calculated system response was initially investigated over a broad 
range of speeds. Results indicated that the critical speeds of the 
system lay between 200 and 900 RPM. Calculations for unit steps of RPM 
were selected so that variations of response with speed would be well 
defined in the critical regions. 

Values initially chosen for the damping constant were somewhat 
arbitrary. Calculations for very small and very large values of the 
damping constant were made. Investigation of system response for all 
possible combinations of damping constants at the two bearings was obvi- 
ously impractical. The case in which the two constants would be kept 
equal seemed to be a suitable starting point. Hereafter the term 
"symmetrical damping constant” or merely "c" will be used to describe 
the stated relationship. 

From the preliminary computer calculations it was noted that the 
most profitable region for close investigation was defined by small 
values of c, i.e., between 0.0 and 3.0 lbs/ (in . /sec . ) . 
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Computer runs were made for no damping, (bearings supported by the 
springs alone), and for damping values up to 3.0 lbs . / (in. /sec . ) in steps 
of 0.5 lbs . / (in . /sec . ) . The case, supposedly analogous to the experimental 
test for damping of 1.77 lbs . (in. /sec . ) , was also run on the computer. A 
value of exceptionally large damping, 10,000 lbs . /(in . /sec . ) , was assigned 
to produce calculated response corresponding to the fixed-bearing exper- 
imental test. Results of the computer calculations described are shown 
graphically in Figs. 24-41. The three modes of vibration associated 
with no damping are shown in Fig. 24. 

4.2.2 Critical and Optimum Damping 

Ordinarily the concept of critical damping is reserved for systems 
with few degrees of freedom in free vibration. In the simple, single 
degree of freedom, spring-mass- da shpot system as treated for example by 
Burton, ^12], the differential equation of motion is developed as follows: 

nmX = - c X — x assuming X = a 3 £ -h ^ - o 

/rn /Til 

c - ~ C ^ ~ 4 nm -k, 

InrTV 



For c -o , 

S~ il 

X = R 

The finite value of c for which the term in the radical reduces to zero 
is referred to as the critical damping constant. From the resulting 
expression for the response of the system at critical damping it can be 
seen that the system does not oscillate; any initial displacement from 
equilibrium is followed by a direct return of the mass to the equilibrium 
position. Larger values of damping result in a similar response but the 
time to return to the equilibrium position is increased. For infinite 
damping the time for return is infinite, or, in other words, the mass 
never returns to its initial position. The pertinent part of the defin- 
ition of critical damping is that it specifies the minimum value of 
damping for which the mass does not oscillate in response to a disturbance. 
For physical systems damping is always greater than or equal to zero. 



POP C£/r/CAL DAMP/ A/ 6 } 

c z =-4/mA c = 2. 

C = — £- 

d 2.001 

__C £ 

x= Re 



- 53 - 



The single degree of freedom system is described by a quadratic 
in 6. If this treatment is extended to a system with two degrees of 
freedom with both masses subject to damping a quartic expression in 5 
results. In this case two damping constants govern the nature of the 
imaginary parts of the two pairs of complex roots. Generally then, for 
a system of n masses there will be an expression of degree 2n in 6 and 
n damping constants will govern the nature of the n conjugate pairs of 
complex roots . 

Applying this treatment to the system under investigation results 
in a sixth degree expression in &. 

/'>/'$ Mr- yU, ->■ ^ = ^ V=co, z -f-oji/4 co, z =ojJ- Z = col 

+ Sr* S** +K)l 4 ‘+(2V+i / L T )& 3 +(V z +Vl/t T +ZX--Jsr --it ) £ Z 

+ (Zvz-Z. 2 -) S + v 2 £ + lt -\Lit =o 

Z !<o 2 

Certain properties of the roots can be predicted from physical consid- 
erations alone. Since only two of the masses are subject to damping 
there will be only two governing damping constants. One mode of oscil- 
lation will be present no matter what values of damping are assigned. 

For certain values of damping the imaginary parts of two of the three 
complex pairs of roots will vanish. 

An extensive analysis of the sixth degree expression is beyond the 
scope of this study. The computer was used in a limited analysis. Results 
indicated that, for the masses and stiffnesses involved in the present 
tests, the imaginary parts of two of the pairs of roots would vanish for 
damping constants of about 30 and 5.0 lbs . / ( in . /sec . ) respectively. For 
these values the system, in free vibration, would exhibit a response some- 
what analogous to a less complex system with critical damping. 

Figures 24-41 show the calculated response of the system over a range 
of small damping constants and for one extremely large value. 
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Figure 24 shows the modes of the three critical speeds for the 
undamped case. The calculated amplitudes of the masses at these three 
frequencies are not infinite because the exact values at which they occur 
were not assigned in the speeds selected for calculation. 

As c is increased from 0.0 to 2.5 lbs . / (in . /sec . ) the amplitude response 
of the masses at the first critical speed is seen to decrease. See 
Fig. 41. Responses at speeds higher than the first critical speed become 
less and less pronounced as damping is increased Increasing the damping 
has a relatively greater effect on the amplitudes at critical speeds 
higher than the fundamental; once the damping exceeds 2.5 lbs . / (in . /sec . ) 
there are virtually no critical speeds higher than the fundamental 

Comparison of Figures 35 and 37 shows a reversal in one trend of the 
response for increasing damping. Although the amplitude response of the 
bearings continues to decrease with higher damping, the rotor response 
is larger at 3.0 than at 2.5 lbs . / (in . /sec . ) . For the strongly damped 
case shown in Fig. 39 the rotor response has increased greatly over that 
for the lightly damped cases. 

Figure 41 shows the amplitude response of the masses as a function 
of the damping constant. It can be seen that the value of 2.5 lbs./ 

(in. /sec.) constitutes that value of c for which the amplitude response 
of the system is a minimum. From this it is apparent that, under these 
conditions, maximum energy is being dissipated from the system. This 
case is defined as that of optimum damping. 

4.2.3 Discussion of Theoretical Results 

Figures 40 and 41 show the major effects of damping on the system 
Figure 40 shows that the critical speed depends directly on the magnitude 
of damping. Similarly, Fig. 41 shows the dependence of amplitude response 
on damping. 

The limiting case of damping on critical speed is that for which 
damping is infinite, or, in other words, fixed supports. The system is 
then reduced to one of the single mass of the rotor. The calculated value 
of this critical speed is 392 RPM. 



- 73 - 



In order to complete the description of the system at each value 
of damping, phase angle diagrams are included with each figure showing 
amplitude response. Somewhat contrary to the behavior of the spring- 
mass-dashpot system, the system considered shows a marked dependence of 
critical speed on damping, and the amplitude response is not always in 
inverse proportion to the damping. 

4. 3 Comments on the Theoretical Analysis 

In general, the theoretical results appear to be valid and consistent 
with established theory. In confirmation of this, a simple, if limited, 
check on the method of calculation is presented. 

For the case of fixed supports, Prohl , f 2j, 

= first critical speed, RPM 
k = shaft spring constant, lbs. /in. 

W = rotor weight, lbs. 

= 187.7 (k/W)^ RPM 
= 187.7 (66.5/15.21)^ = 392.3 RPM 
The calculated critical speed for this configuration was 392 RPM. 
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5. COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 



Figure 31 shows the calculated and test amplitude responses for the 
case of damping of 1.77 lbs . / ( in . /sec . ) . 

The observed critical speed was 383 RPM From Fig. 40, First 
Critical Speed vs. c, the critical speed corresponding to c 3 1.77 lbs./ 
(in. /sec.) is 346, a discrepancy of 37 RPM. Also from Fig. 40, a critical 
speed of 383 RPM corresponds to a damping factor of 6.5 lbs. /(in. /s ec.). 
From analysis of the observed critical speed there is the indication that 
the damping during the test was greater than that of 1.77 lbs . / ( in . /sec . ) 
thought to have been in effect. 

Amplitude indications from the test were analyzed in conjunction with 
Fig. 41, Mass Amplitude at First Critical Speed vs. c. From the figure 
it can be seen that an amplitude of 0-116 corresponds with either of 
two damping constants, 0 047 or 18.0 lbs . / (in . /sec . ) . The lower value is 
felt to be totally improbable. The higher value, although excessive in 
itself, tends to confirm the indication noted in comparison of critical 
speeds that the damping was larger than 1.77 lbs . / (in. /sec . ) . Amplitude 
measurement error, although possibly greater at critical speeds due to 
transients from non-steady drive, could not account for the large dis- 
crepancy . 

The apparent discrepancy in dashpot damping could stem from any of 
several sources. The dashpots may have been cocked during the run 
and produced a strong non-linear response. The dashpots were mounted 
vertically for their calibration. Their operation in the horizontal 
plane during the critical speed test may have produced the apparent 
change in their characteristics. 

Figure 39 shows the calculated and experimental response for the 
fixed-bearing case. The calculated critical speed was 392 while the 
experimental value was 380. Since amplitudes are theoretically infinite 
for the fixed-bearing case, comparison is not applicable as in the case 
with lightly damped bearings. Based solely on the comparison of critical 
speeds the agreement between calculated and test results is fair. 
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6 CONCLUSIONS 



The system of a rotor mounted on a shaft with damped and flexibly 
supported end bearings shows its amplitude response to be dependent on 
bearing support impedance. Rotor amplitudes decreased as damping in- 
creased to the optimum value of 2.5 lbs . / (in . /sec . ) and increased for 
larger values of damping. There is a limit to which rotor amplitude 
at critical speeds can be influenced by damping at the bearings and one 
critical speed will persist in spite of damping applied. 

The calibrated damping characteristics of the dashpots were either 
in error or dashpot damping was effected by the mechanical arrangement 
of the dashpots on the critical speed rig. 

Accuracy of test measurements could be improved. 
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7. RECOMMENDATIONS 



The critical speed test rig offers promise for use in further 
vibration analyses if improvements in its damping and control charac- 
teristics can be effected. 

It is recommended that the physical system be improved by 

a. Modification of the drive system to eliminate thrust along the shaft 
axis , 

b. Design of a support for spring-spring holder assembly to permit testing 
without dashpots, 

c. Provision of springs with stiffnesses between those of the springs 
tested, and 

d. Use of better speed and amplitude indicators. 

It is also recommended that a further investigation of the damping 
range and characteristics of the dashpots be conducted. 
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